Structural and magnetic dynamics of a laser induced phase transition in FeRh 
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We use time-resolved x-ray diflfraction and magnetic optical Kerr effect to study the laser induced 
antiferromagnetic to ferromagnetic phase transition in FeRh. The structural response is given by 
the nucleation of independent ferromagnetic domains (ri ~ 30ps). This is significantly faster than 
the magnetic response {t2 ~ 60ps) given by the subsequent domain realignment. X-ray diffraction 
shows that the two phases co-exist on short time-scales and that the phase transition is limited by 
the speed of sound. A nucleation model describing both the structural and magnetic dynamics is 
presented. 
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To date the fastest manipulation of magnetic films and 
elements are induced by single fs laser pulses and include 
domain switching and demagnetization on sub-ps time- 
scales [ll-Hl- On similar time scales magnetization has 
been switched by ultrashort but strong magnetic field 
pulses generated by relativistic electron bunches and 
stripline techniques 0, while a third intriguing option 
is the manipulation of the magnetic energy landscape 
by strong single cycle electric field pulses [Tf. Gener- 
ation of a magnetic moment is equally interesting but 
harder to achieve on an ultrafast timescale. Ferromag- 
netic (FM) order can be established in ordinary FM ma- 
terials by cooling from the paramagnetic phase, but the 
process is limited by heat transfer and typical timescales 
are nanoseconds. In this context the antiferromagnetic 
(AFM) to FM phase transition in FeRh (Tt « 375K) is 
interesting because it can be induced by fs laser pulses. 
The first order phase transition from the low tempera- 
ture AFM phase is accompanied by a '^0.5 % volume 
increase, and though this has been known since 1939 [g] 
the physical mechanism behind the transition has never 
been resolved and is still debated 

The ultrafast transition has been studied in all-optical 
pump-probe experiments using both transient reflectiv- 
ity and time resolved magneto optical Kerr effect (TR- 
MOKE) [H [13 . While the reflectivity measures a 
combination of electronic and structural properties, TR- 
MOKE measures the magnetization. In addition x-ray 
magnetic circular dichroism (XMCD), which is element 
specific, has been used to probe the transition 14|. This 



study found a gradual growth of the magnetization on a 
time-scale of ~ 100 ps. In none of these experiments was 
a separate determination of the lattice dynamics possi- 
ble. The magnetization dynamics have been simulated 
with the Landau-Lifshitz-Gilbert equation and a model 
which considers growth of a non-homogenous magneti- 
zation proportional to the spin-temperature followed by 
a slower alignment and precession of the local magnetic 



moments |15l |. 

In order to unravel the physical processes underlying 
the phase transition a prerequisite is the ability to dis- 
tinguish the contributions arising from the lattice- and 
magnetic- dynamics. In this letter we report a time re- 
solved x-ray diffraction (XRD) experiment directly mea- 
suring the structural dynamics and providing the evolu- 
tion of the AFM and FM volume fractions with a time 
resolution of ^ 200 fs. The results are combined with 
TR-MOKE measurements on the same sample and are 
explained using a simple model describing the nucleation 
and subsequent alignment of FM domains. 

The FeRh epitaxial thin film (d = 47 nm) was grown 
on MgO(OOl) by co-magnetron sputtering from elemen- 
tal targets. The film is epitaxial with a (001) surface. 
Upon heating through the phase transition the lattice ex- 
pands 0.7 % along the surface normal, in contrast to the 
isotropic expansion of bulk samples, indicating a strong 
in-plane strain. 

The time-resolved XRD was performed at an x-ray 
energy of 7 keV using a synchrotron slicing source 
(200 ph/pulse at 2 kHz and 1.2 %bw) as probe and an 
800 nm 120 fs p-polarized laser pulse with an incidence 
angle of 12° as pump, resulting in a total time resolu- 
tion of 200 fs 16]. The x-ray gracing incidence angle 
was either a = 0.51° or 0.71°, in order to match the 
penetration depth of the x-ray probe to either the laser 
pump (penetration depth 15 nm) or to the film thick- 
ness (47 nm). Due to the grazing angle the x-ray spot 
size was 0.4 x 1 mm^. The (101) Bragg refiection was 
recorded with a two dimensional PILATUS lOOK pixel 
detector and rocking curves with 60 discrete images were 
acquired by rotating the sample around the surface nor- 
mal (±2.5°) [13. 

The TR-MOKE measurements were performed at 
72 kHz as a two-color pump-probe experiment with 200 fs 
cross-polarized 800 nm p-polarized pump and 400 nm 
probe pulses. The probe (58 x 26 /im^) had an incidence 
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FIG. 1: (Color online) Shift of the center of the (101) Bragg 
peak after laser excitation. (A) As a function of time delay. 
(B) As a function of fluence at a fixed time delay of 145 ps 
for temperatures above and below Tt ~ 375K. 



angle of 60° with respect to the surface while the pump 
was incident along the surface normal (90°). We used a 
longitudinal geometry with an applied in-plane external 
field of 0.1 T. The time dependent traces were recorded 
for opposite orientations of the external field _and the Kerr 
rotation is the difference of the two traces 
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In FIG.[B\ we show the shift of the center of the Bragg 
peak as a function of time and laser fluence. As the ex- 
pansion due to the film geometry is purely one dimen- 
sional, the shift occurs along the surface normal {q±). 
A negative peak shift corresponds to an expansion of the 
lattice, and for a single unstrained phase the peak shift is 
proportional to the lattice expansion, Aa/a w —Aq±/q±. 
We first consider data obtained at a sample temperature 
of T = 340 K < Tt. At a pump fluence of 1 mJ/cm^ the 
Bragg peak shifts to a new equilibrium position through 
a single damped oscillation. This is the signature of a 
thermally induced strain wave [l9| and the period T— 
18.6±0.9 ps (95% confidence interval) is given by the time 
it takes a strain wave to travel back and forth through 
the film, T = 2d /v. The resulting speed of sound v — 
5.1 km/s is consistent with the literature jl^]. At higher 
pump fluences tlie strain wave is still present but ac- 
companied by an increased peak shift. At intermedi- 
ate (2-2.7 mJ/cm^) laser fluences this extra peak shift 
is slower than the strain wave, while at high laser flu- 
ences (5.2 mJ/cm^) the time scales are comparable. We 
attribute this peak shift to the transformation from the 
AFM to the FM phase, and since the peak shift at in- 
termediate fluences is slower than the strain wave, the 
strain wave can not be the driving force of the transi- 
tion. In FIG. [Tj3 the peak shift is shown as a function of 
fluence for a time delay of 145 ps. A deviation from the 
linear fluence dependence given by thermal expansion, 
is observed above 1 mJ/cm^. This threshold behavior 
is the signature of a laser-induced phase transition. The 
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FIG. 2: (Color online) Bragg reflections obtained at different 
time delays, (A) for a pump fluence of 5.2 mj/cm^ and (B) for 
a pump fluence of 2.7 mj/cm^. The solid line is the simulated 
Bragg reflection at t=6ps under the sole influence of thermal 
strain. 



threshold fluence of ~1 mJ/cm^ corresponds to a temper- 
ature increase of 32 K [2^. This matches the difference 
between the sample (340 K) and the transition temper- 
ature (Tt=375 K) and is consistent with the thermal 
nature of the laser induced transition. As a final verifi- 
cation that the peak shift arises from the AFM to FM 
phase transition the sample was heated to 440 K, well 
above T^. At this temperature the peak shift is given 
solely by thermal expansion. We thus unambiguously 
confirm the laser induced AFM to FM phase transition. 

In FIG. [5J\ we further show the full Bragg peaks 
for several time delays and a fixed pump fluence of 
5.2 mJ/cm^. At i = ps the entire film is in the AFM 
phase. At t — 400 ps the entire film has been driven into 
the FM phase. The full transformation is evident as the 
two peaks have the same shape. At intermediate times 
the peak is a sum of the two distinct peaks corresponding 
to the two phases. We thus directly confirm the coexis- 
tence of the two structural phases at short time scales. 
This is supported by the solid line which is the calcu- 
lated Bragg peak at t = 6 ps under the sole influence of 
thermal strain J^] . It is thus evident that the transition 
proceeds through the decrease of the AFM phase and in- 
crease of the FM phase, rather than through a continuous 
change of the lattice constant. The phase coexistence is 
a common trait of first order transitions and has been 



observed statically in FeRh for both the magnetic 
and the crystallographic structure 



21|. FIG. [2^ shows 



similar Bragg peaks obtained at a fixed pump fluence 
of 2.7 mJ/cm^. In this case the final state {t = 400ps) 
consists of both AFM and FM structure. 

In FIG. [IK it can be seen that the amplitude of the 
initial strain wave does not scale with fluence above the 
threshold. Thus the phase transition adds to the mag- 
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nitude of the strain wave. The fact that the phase of 
the FM peak is initially strained is also seen in FIG. [2K. 
Since the stress is relaxed at the speed of sound at least 
for parts of the film the underlying change in energy land- 
scape from AFM to FM appears faster than r ^ d/v. 
This timescale is the limit for the structural transition. 

To obtain a quantitative measure of the structural 
change due to the phase transition we separate it from 
the effects of the strain wave. In our experiment the 
grazing incidence geometry and the relatively low x-ray 
energy limit the separation Agj_ between the AFM and 
FM peaks, but the two contributions can be systemati- 
cally extracted from the measured data by fitting it to 
the sum of two symmetric functions (a/ cosh((x — b)/cY) 
[26| . One such fit is shown in the insert in FIG. 
This way the integrated intensities which are propor- 
tional to the scattering volumes of the two phases can 
be extracted, and in FIG. [3K the volume fraction of the 
FM phase (W fm) is shown as a function of time delay. 
Before the arrival of the laser pump the film is entirely in 
the AFM phase and Vfj\/=0. After laser excitation the 
FM volume grows to saturation within 100 - 200 ps. In- 
creasing the fluence from 2.0 to 2.7mJ/cm^ increases the 
transformed volume fraction of the FM phase ^*pM from 
32% to 59%. By decreasing the x-ray penetration depth 
from 47 nm (a = 0.71°) to 15 nm (0.51°) we see that 
for smaller probe depths Vfm rises significantly faster, 
while the transformed volume is only increased slightly. 
This implies that the nucleation of the FM phase starts 
at the free surface while the final phase consists mainly 
of domains penetrating the entire film depth. 

To compare the change in structure to the change in 
magnetization we measured the change in Kerr rotation 
(AOfcerr) On the same sample but at 313 K, as shown in 
FIG. 03. The higher pump fluences applied compensate 
for the lower sample temperature. As for the x-ray data 
we observe a threshold in laser fluence. Below the thresh- 
old there is no change in magnetic moment (1.3 mJ/cm^) 
while above the threshold the Kerr rotation increases 
and reaches saturation after several hundred ps (3.1 and 
4.5 mJ/cm^). For the first 100 ps the dynamics appears 
to follow a power law, in strong contrast to the growth of 
Y FM- As shown in the model below this difference arises 
because XRD measures the volume of FM domains, while 
TR-MOKE measures their alignment. Our TR-MOKE 



results are in agreement with previous XMCD [ij] and 



TR-MOKE [1^ [13[ experiments, except for the previ- 



ously reported ultrafast TR-MOKE component which we 
do not observe. The absence of the ultrafast response 
may have two reasons. First, TR-MOKE is prone to op- 
tical artifacts in the first hundreds of fs after excitation. 
Second, the absence might be due to the different probe 
spot size used in the different experiments, as this deter- 
mines whether or not the magnetization is averaged over 
many FM domains or just a few. 

We now present a model describing the observed dy- 
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FIG. 3: (Color online) (A) Evolution of the volume fraction 
of the expanded FM phase as a function of time delay, shown 
for two different pump fluences and incidence angles. The 
insert illustrates the fitting procedure (t = 6ps, 2.7 mj/cm^), 
showing the AFM (dashed) and FM (dotted) fit functions, and 
their sum (solid). (B) Transient Kerr rotation for comparable 
fluences. The inset shows the dependence of the time constant 
ri (o) and a chi-square estimator (+) on the model parameter 
n. The solid lines are fits to the data as described in the text. 



namics in terms of nucleation and alignment of FM do- 
mains. We assume that the film is instantaneously heated 
above Tt and that the nucleation of the FM phase pro- 
ceeds through nucleation at many independent sites. The 
rate of change of Vfm is then proportional to Vafm and 
described by a single time constant ri which may differ 
for the structural and magnetic changes. Because the 
transition is of first order the final state may be mixed. 
To account for this we introduce the final fraction of FM 
phase YpM- ^ fm + ^ afm = 1 we then find: 



FM 



Vafm - (1 - V^j 



FM) 



dt 



Tl 



^FAf ~ ^FM 
Tl 



(1) 



This is essentially the Avrami model without growth of 
existing domains (2^ and the solution is the exponential 
function V^a/(1 — exp(— t/n)). The depth dependence is 
only included by allowing different time constants when 
averaging across different probe depths. This exponential 
function has been used to fit the data in FIG.[3l4. and de- 
scribes the data very well. Within error bars we find the 
same time constant ri = 33 ± 4 ps for both intermediate 
fluences, while ri = 14±3psata = 0.51°. The exponen- 
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tial growth of Vfm is only observed when the nucleation 
of independent domains is dominant while the growth of 
existing domains is suppressed. As the x-ray spot size is 
0.4 X 1 mm^ the distance between nucleation sites must 
then be less than ~ 10 ^m. 

In order to describe the magnetization we assume that 
the FM domains nucleate with the magnetization fixed 
into one of n directions and that all n directions are 
equally probable. We assume one of these directions is 
favored by the applied magnetic field and define the vol- 
ume fraction Va of FM phase aligned to the applied field 
and the volume fraction V„ not yet aligned. These sat- 
isfy Vfm = Vn + Va- We finally assume that alignment 
of the magnetization occurs by growth of the aligned FM 
domains at the expense of non-aligned FM domains, as 
described by a product term V^V„/t2 and a single time- 
constant T2. As this model depends on the existence of 
both aligned and un-aligned FM domains it supports a 
theory where short-range interactions are responsible for 
aligning neighboring FM domains through domain wall 
motion. Given these assumptions the time evolutions of 
Vn and Va are described by two differential equations: 

dVn _ n-1 dVpM _ VaVu 
dt n dt T2 

dVA _ 1 dVpM ^ VaVu 
dt n dt Ti 

The MOKE signal is proportional to the magnetization 
along the preferred direction: <m> cx - V„/(n-l). 
This is valid when one of the n-1 directions is opposite 
to the applied field and the magnetization of the remain- 
ing n-1 directions averages to zero. The underlying as- 
sumptions of nucleation and coexistence of phases in this 
model differ from the work by Bergman et al. [l^ who as- 
sumed that the local magnetization grows monotonously 
with spin-temperature throughout the film. 

The three differential equations are solved numerically 
for integer values of n, and the result of fitting the result 
to the TR-MOKE data is shown in FIG. [3^ for 7i=4, 
which optimizes the fit. The agreement between experi- 
ment and fit is excellent with = 1.8. The result n=4 
is consistent with the in-plane magnetization expected 
for a cubic thin film, and has been confirmed by static 
XMCD PEEM images obtained for the same film. For 
the nucleation time t\ = 17.8 ± 0.9 ps we find, as for the 
XRD data, the same value for both pump fluences within 
errorbars. Since the probe depth for the laser in the TR- 
MOKE experiment is '^ll nm this must be compared to 
the 14 ± 3 ps obtained with XRD at a = 0.51°. We thus 
conclude that the timescales of nucleation for magnetic 
and structural domains are the same within the errorbars. 
For the final parameter T2 which describes the alignment 
of domains we obtain 72 ± 1 ps and 57 ± 1 ps for fluences 
of 4.5 mJ/cm^ and 3.1 mJ/cm^ respectively. Based on 



the good agreement between data and model we conclude 
that the FM domains initially nucleate with un-aligned 
moments which are subsequently aligned. We speculate 
that the initial domain structure is given by the underly- 
ing AFM phase and that the mechanism responsible for 
the re-alignment process is domain wall motion. 

The simplest alternative model would describe realign- 
ment as a rotation of the total moment of a domain. The 
realignment term in ([2]) and ([3]) would then be TV,i/r2, 
independent of Ya- For this model w 30 which is 
significantly worse than = 1-8 obtained above. In ad- 
dition the fit is independent of n and results in different 
nucleation times ti for structure and magnetism. This 
would imply that the phase transition is independent of 
the anisotropy and that the magnetic domains nucleate 
slower than the structure. Both conclusions appear less 
likely than those reached from ([T]) - We thus reject 
the alternative explanation. 

In summary we have measured the structural and mag- 
netization dynamics of the AFM to FM phase transition 
in FeRh on an ultra-fast timescale. XRD allowed us to 
directly observe the co-existence of the two phases and 
to derive a simple model which describes the evolution of 
both the structure and the magnetization. We find two 
intrinsic timescales: One for the initial nucleation of FM 
domains which is the same for both magnetic and struc- 
tural dynamics, and a second for the subsequent growth 
of FM domains aligned to the applied magnetic field. The 
co-existence of FM and AFM domains is clearly seen in 
the XRD data. At intermediate pump fiuences the phase 
transition to a large extent proceeds similarly to static 
heating, while at higher fluences the structural change 
is limited by the speed of sound. This speed limit on 
structural change in principal allows for a significantly 
faster magnetic response, which we do not observe. It 
thus appears that magnetic and structural nucleation go 
hand in hand rather than one driving the other. While 
the microscopic nature of the magnetization change has 
been considered theoretically [lH, a more definitive an- 
swer will require the use of spatially resolved magnetic 
probes in order not to average the dynamics over many 
domains. 
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